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Scaling functions for Tsallis non— extensive statistics. 
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We study the one-dimensional Ising model with long-range interactions in the context of Tsallis 
non-extensive statistics by computing numerically the number of states with a given energy. We 
find that the internal energy, magnetization, entropy and free energy follow non-trivial scaling laws 
with the number of constituents N and temperature T. Each of the scaling functions for the internal 
energy, the magnetization and the free energy, adopts three different forms corresponding to q > 1, 
q — 1 and q < 1, being q the non-extensivity parameter of Tsallis statistics. 
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It is generally assumed that Thermodynamics and Sta- 
tistical Mechanics necessarily imply that the entropy, the 
internal energy and other thermodynamic potentials are 
extensive quantities. For instance, the internal energy 
E as a function of temperature T and number of con- 
stituents N scales usually as: 



E(N,T) =Ne{T). 



(1) 



Within the theoretical framework of Statistical Mechan- 
ics, this is indeed a widespread consequence when the in- 
termolecular potentials are short-range. For the so called 
normal systems M, the number of microscopic states 
with a given energy scales as £l(E, N) — exp(Ns(E/N)), 
from where it follows the entropy behavior S(E, N) — 
In ft = Ns(E/N) J|. The thermodynamic relation 
T _1 = (f^)jY leads then to the scaling law ([[]) for the 
internal energy. However, it has been also realized that 
long-range potentials can lead to non-extensive behav- 
ior and, recently, there has been some interest in finding 
the correct scaling laws for the thermodynamic potentials 
for systems whose non-extensive behavior arises from a 
long-range interaction [0,0. Let us be more specific and 
consider the ferromagnetic Ising model with long-range 
interactions: 



A ! 



H ~ 1^ —^ 



(S, = ±1, V»), (2) 



i,J=l 



1,3 



where indexes i, j run over the N sites on a d-dimensional 
lattice and r^j is the distance between sites i and j. It 
can be easily shown that the energy levels scale as 



NN* = N- 



1-a/d 



(3) 



In the case of a > d, N* tends to a constant in the limit 
of large N and the energy recovers its usual extensive 
behavior, whereas in the case a < d the behavior is non- 
extensive (for a = d the limit N* = InN is assumed). 
Therefore one expects the failure of the scaling law p) 
for a < d. This is indeed the case as Cannas and Tamarit 



B have shown by performing Monte-Carlo simulations of 
the Hamiltonian (0) in a d=l system. Their results show 
that the Boltzmann-Gibbs canonical ensemble statistics 
leads to the following scaling laws for the internal energy, 
spontaneous magnetization, entropy and free energy: 



E(N,T) = NN*e(T/N*), 
M(N,T) =Nm(T/N*), 
S(N,T) = Ns{T/N*), 
F{N,T) = NN*f{T/N*). 



(4) 
(5) 
(6) 

(7) 



The argument justifying these scaling laws can be sum- 
marized as follows pj: the internal energy and the en- 
tropy appear in the definition of the Helmholtz free en- 
ergy as F = E — TS, therefore one expects that E and 
TS should have the same behavior for large N. Since E 
scales as NN* and S scales as N one obtains that T must 
scale as N* thus leading to the previous scaling ansatzs. 

Although the above scaling laws have been verified in 
H by application of the Boltzmann-Gibbs statistics, it 
has been argued that the appropriate frame to describe 
systems with long-range interactions should be that 
of Tsallis non-extensive statistics, since non-extensivity 
properties appear in this formulation in a natural way J5| . 
Tsallis statistics depends on a parameter q in such a way 
that the limit q = 1 retrieves the results of Boltzmann- 
Gibbs statistics whereas for q < 1, the entropy is super- 
extensive and for q > 1 it is sub-extensive. 

The aim of this paper is to derive and compute numer- 
ically the scaling laws for the entropy, internal energy, 
free energy and magnetization that follow form the ap- 
plication of Tsallis statistics to the long-range d=l Ising 
model defined by (0) in the non-extensive regime a < d. 
Our main result is that we can write scaling laws (see 
Eqs . (p5|)— (p8|) below) that depend on appropriate scal- 
ing factors Ag(N), Af(N), A%(N) and N*. In the limit 
q —> 1 the scaling laws for Boltzmann-Gibbs statistics 
(0)-(n) are recovered. Furthermore, the scaling functions 
e q , m q and f q depend on the parameter q in such a way 
that they collapse onto only three scaling functions for 



each magnitude: those of q > 1, q = 1 and q < 1. 

Let us remind briefly which are the basic ingredients of 
Tsallis statistics. Each of the W system configurations 
(W — 2 N for the Ising model used here) is assigned a 
probability pi , which is obtained by finding the extrema 
of the generalized entropy 



S q = 



1 



-T w v q 



(8) 



subject to appropriate constraints. Once the p^s have 
been obtained, the quantities of interest are computed as 
generalized averages of microscopic functions Oi M : 



(0) g 



L^i=\ Pi 



(9) 



In the canonical ensemble, the main constraint (be- 
sides the normalization condition E».P» — 1) is that the 
mean value of the energy is fixed to a given value E q . 
This variational problem has the implicit solution for the 
configuration probabilities: 



Vt = 



[l-(l-q)/3'e t } — 



*,w 



Ejiji-a-?)^-]* 



i i 



(10) 



where ej is the energy of the «-th configuration. We have 
used the notation 



0' = 







w 



(i - q)P Ei=i phil Ei=i p! + Ej=i p. 



™ *) 



(ii) 



and the Lagrange multiplier = 1/T (the equivalent of 
the inverse temperature for the Boltzmann-Gibbs canon- 
ical ensemble) has to be found by imposing that the 
mean value of the Hamiltonian is equal to the given value 
E q = (Ti.) q . The usual procedure, however, is to give a 
value for T — 1/0 and to derive, using equations Jig ) 
and (O), the probabilities Pi{0) as a function of the (in- 
verse) temperature and then compute the mean value 

EM = Y,ZiPiWei/TZiPi(P) q - 

Now the main problem arises in Tsallis statistics that 

we do not have the solution for the probabilities Pi(0) 
in a closed form, since the non-linear coupled equations 
(|l(j)-([ll]) have no explicit solution. Of course, in the case 
of q = 1 we do know the solution (up to a normalization 
factor) which is nothing but the celebrated Boltzmann 
factor: Pi((3) = Z~ 1 e~" ei where Z is the partition func- 
tion. The explicit knowledge of the probabilities pi (0) in 
the case q — 1 allows the use of Monte-Carlo techniques 
for the numerical calculation of the averages (0). In its 
simplest version H, the Metropolis algorithm proposes 
a new configuration j by randomly flipping one spin in 
configuration i. The new configuration j is accepted with 
a probability vam.(l,Pj/pi) — min(l,e _/3 ^ j_£i ' ) ). Notice 
that the partition function cancels out in the calculation 



of the acceptance probabilities. Unfortunately, since for 
q =/= 1 the probabilities pi are not known as a function of 
0, there is no trivial generalization of the Monte-Carlo 
method to perform the averages in (O) at fixed temper- 
ature 0. One can perform Monte-Carlo simulations at 
fixed 0' M, but then the physical temperature is not 
known. Another interesting approach (close in spirit to 
our method here) is that of Lima et al. M who have used 
the broad histogram method Il0[ to study the 2-d short- 
range Ising model, focusing mainly on the possibility of 
the existence of a phase transition for q =/= 1 . 

We overcome these problems by using a method of 
histogram by overlapping windows initially devised to 
study short-range lattice models [[llj. In this method, 
one computes numerically the number Q(E, N; SE) |1J] 
of microscopic states whose energy lies in the interval 
{E,E + SE). The histogram by overlapping windows 
method performs a microcanonical simulation by fixing 
the energy in a window {E, E + AE) and computing 
the ratios n(E u N; SE)/n(E 2 ,N; SE) for energies E t , E 2 
within this window. Once those ratios have been com- 
puted with a given accuracy, we perform another micro- 
canonical simulation in a different window (£", E' + AE) 
which overlaps the previous energy window. The method 
proceeds until the windows have swept over all the possi- 
ble energy values. The exact knowledge of the degeneracy 
for the ground state £l(E ,N) — 2 allows the recursive 
calculation of the number of states il(E, N; 5E) for all 
values of E. For the long-range Ising model, the size of 
the window AE has to been chosen carefully in order to 
avoid the lack of ergodicity. A full account of the method 
details will be given elsewhere [[L3|. Here we just report 
on the results we obtain for the aforementioned scaling 
laws. 

Using this method we have computed the number of 
states Q(-Efc) for the d=\ Hamiltonian defined in (0) with 
a=0.8 and system sizes N = 34, 100, 200, 400, 1000. 
Once the number of states ft(Ek) is known, one can use 
a recursive method p5l to solve Eqs.([J7)|)-(|ll"|) in order to 
find the probabilities pi (/3) . Equivalently, one can com- 
pute the probabilities Pi(/3') as a function of the param- 
eter /3' using (|10|), where the sum over configurations is 
now replaced by a sum over all possible energy bins of 
size SE. The entropy, S q (f3'), the internal energy E q (j3') 
and the magnetization M q ((3') [[14| are computed in the 
same way as a function of (3' using relations (g) and (|j). 
The physical temperature T = 1/(3 can be obtained by 
inverting ([il|): 



„o, 1 - (g - W 

' ' l-(l-a)0'E„(0'Y 



(12) 



thus allowing a parametric plot of the internal energy, 
entropy and magnetization as a function of the temper- 
ature T. We have also computed the free energy defined 
as F q = E q — TS q . It is important to remark that, in the 



case of q < 1, and for large values of system size N the 
raw data show a loop with temperature. This is similar 
to what happens in the short-range Ising model and we 
have adopted the same criterion than in Jig] : to use a 
Maxwell-type construction that replaces the loop of the 
curve by a straight line joining the two points with the 
same value of the free energy. 




Keeping in mind that the energy scales as NN* and as- 
suming that E q and TS q scale in the same way we de- 
rive that the temperature must scale as NN* /A q (N) and 
hence we are led to the ansatz: 



E q (N,T) = NN*e q (TA^(N)/NN*), (15) 

M q (N,T) = Nm q (TAf(N)/NN*), (16) 

S q (N, T) = A q (N)s q (TA s q (N)/NN*), (17) 

F q {N,T) = NN*f q (TA q (N)/NN*). (18) 
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FIG. 1. Number of states Q(E,N) plotted to check the 
scaling law ([f3J). The results for N — 34 have been obtained 
by an exact enumeration of the W — 2 34 possible states, 
whereas the results for the other system sizes have been ob- 
tained by the histo gram by overlapping windows method de- 
scribed in the text J16J. 
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FIG. 2. Internal energy plotted to check the scaling relation 
(151) by using different values of q and system sizes N. The 
scaling factors used are those defined in (hj) and (M). 



Here, in view of later results, we have introduced new 
scaling factors A q (N) and A q (N). The previous argu- 
ment would imply simply Af(N) = A q (N) = A q (N). 
Notice that in the limit q — ► 1 it is Ai(N) ~ TV and the 
scaling laws Eqs.(Q)-(0) are recovered. 

In figures (pi) and (0) we scale the energy, magnetiza- 
tion and entropy data by using factors A q (N), A q (N), 
A$(N) and N* as implied by Eqs.dH)-©. In figure 
(g) we concentrate in the validity of scaling for different 
values of N, whereas in figure (g) we compare the scaling 
functions for different values of q using the scaling func- 
tions obtained for the largest value N — 1000. These 
figures give evidence that in the case q < 1, scaling is 
well satisfied by using Af(N) = A%(N) = A q (N) as ar- 
gued before. However, the data for q > 1 do not follow 
this scaling description. In order to obtain a good scaling 
for q > 1 one observes numerically that it is necessary to 
assume instead the limits A^(N) ~ 2 N{1 - q ^/{q - 1) and 
A%(N) - 2 N( -i- i y(q - 1). A unifying description that 
reproduces the required limits for all values of q is: 



A s q (N) 



2 N\l-q\ _ 1 

ll-al : 



Af(N) 



A q {N? 
A$(N) 



(19) 



In Fig.(|jj) we plot the computed values for the number 
of states il(E, N) as a function of the energy E for dif- 
ferent systems sizes. In this figure, the data have been 
scaled to show that the number of states follows the scal- 
ing law that one would expect given that the total num- 
ber of states scales as 2^ and the energy levels scale as 
NN*, namely: 



n(E,N) = exp(N<f>(E/NN*)). 



(13) 



In order to generalize the scaling functions for the ther- 
modynamic potentials in the case of Tsallis statistics, we 
notice that, in the case of equiprobability (corresponding 
to very high temperature) (@) implies that the entropy 
scales as S q (N) ~ A q (N), where 



A q (N) = 



1 _ 2 JV(i- 9 ) 

q~l 



and these expressions have been used to scale data as 
shown in the figures. We observe, see Fig.(g) for the inter- 
nal energy, that the quality of the scaling is rather good 
and improves, as expected, with increasing system size. 
A very interesting feature is that, as shown in Fig.(gj) the 
scaling functions group into three different forms corre- 
sponding to q < 1, q = 1 and q > 1. The only exception 
is that of the entropy for which the collapse for q > 1 is 
very poor. This is easily understood by noticing that the 
low temperature limit of the entropy for infinite system 
size is S q (T = 0) = (1 - 2 1 -i)/(q - 1) whereas the high 
temperature limit is S q (T — > oo) = l/(q — 1) and those 
two finite values can not be rescaled simultaneously. The 
scaling for the free energy follows directly from its defini- 
tion F q — E q —TS q . For q < litis/ g (x) = e q (x)—xs q (x), 
(14) whereas for q > 1 and in the limit of large N, the scaling 



function is given simply by f g (x) = e q (0) — xs q (oo) 
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FIG. 3. Internal energy (top graph), magnetization (mid- 
dle graph) and entropy (lower graph) plotted in order to check 
the proposed scaling behavior. We have used in all the curves 
the value N = 1000 and varied the parameter q. For clar- 
ity, in the entropy curve, the insert shows all the values of 
q, whereas the main plot takes only q > 1. The curves with 
q > 1 include q = 1.2, 1.4, 1.6, 1.8 and the curves with q < 1 
include q = 0.2, 0.4, 0.6, 0.8, although the different curves 
are almost indistinguishable with the resolution of this figure. 



In summary, the scaling laws given by Eqs.(|T5|)-([Tq) 
work for all values of q when one uses the scaling factors 
given by Eqs.(Oh, and (fl9|). Moreover, the scaling func- 



tions e 



</• 



and f q adopt only three different forms for 
each magnitude: one valid for q > 1, one valid for q = 1 
and another valid for q < 1. 

Several final comments are in order. First, it is dis- 
tressing the fact that the scaling forms for q > 1 do not 
follow the scaling ansatz that follows naively from the 
argument that T should scale as E q /S q ~ NN*/A q (N). 
We have not been able to find a convincing argument 
that reproduces the scaling forms found in this paper for 
q > 1. It seems that the transformation /?' — ► j3 given by 
Eq.(|l2|) has two special points where the slope changes 
abruptly and which scale precisely as NN* /A^(N) and 
NN*/Aq(N), although the exact implication for the 



scaling functions is not clear to us at this moment. Sec- 
ond, the fact that the scaling functions adopt very dif- 
ferent forms (for q < 1, q = 1 and q > 1) might allow 
to conclude easily whether classical Boltzmann-Gibbs or 
Tsallis statistics should be used when analyzing experi- 
mental data. Finally, we would like to stress the power of 
the histogram by overlapping windows method to study 
numerically systems with long-range forces. 
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